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ABSTRACT 

In 2012 Markarian 421 underwent the largest flare ever observed in this blazar at radio fre¬ 
quencies. In the present study, we start exploring this unique event and compare it to a less 
extreme event in 2013. We use 15 GHz radio data obtained with the Owens Valley Radio Ob¬ 
servatory 40-m telescope, 95 GHz millimeter data from the Combined Array for Research in 
Millimeter-Wave Astronomy, and GeV y-ray data from the Fermi Gamma-ray Space Tele¬ 
scope. The radio light curves during the flaring periods in 2012 and 2013 have very different 
appearances, both in shape and peak flux density. Assuming that the radio and y-ray flares 
are physically connected, we attempt to model the most prominent sub-flares of the 2012 
and 2013 activity periods by using the simplest possible theoretical framework. We first fit a 
one-zone synchrotron self-Compton (SSC) model to the less extreme 2013 flare and estimate 
parameters describing the emission region. We then model the major y-ray and radio flares 
of 2012 using the same framework. The 2012 y-ray flare shows two distinct spikes of similar 
amplitude, so we examine scenarios associating the radio flare with each spike in turn. In the 
first scenario, we cannot explain the sharp radio flare with a simple SSC model, but we can 
accommodate this by adding plausible time variations to the Doppler beaming factor. In the 
second scenario, a varying Doppler factor is not needed, but the SSC model parameters re¬ 
quire fine tuning. Both alternatives indicate that the sharp radio flare, if physically connected 
to the preceding y-ray flares, can be reproduced only for a very specific choice of parameters. 

Key words: galaxies; active - galaxies: jets - BL Lacertae objects: individual: Markarian 
421 


1 INTRODUCTION 

Blazars are active galactic nuclei (AGN) that are variable over the 
entire electromagnetic spectrum. Their variability is enhanced due 
to Doppler-beamed emission from a relativistic jet pointing close 
to the line of sight. Blazars having featureless optical spectra or 
showing emission line s with equivalent w idth < 5 A are historically 
classified as BL Lacs dStocke et al.lll99ll) . 

Markarian 421 (hereafter Mrk 421) is one of the best stud- 
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ied BL Lac objects. It is rela tively nearby with a redshift of 
z = 0.031 jUlrich et alj 1 197 51). It is classified as a high syn¬ 
chrotron peaked blazar lAbdo et alj|2010bl) based on its spectral 
energy di stribution (SEP). It was the first blazar detected at TeV 
energies jPunch et alJll992h . and since then it has been the sub- 
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Tosti et alT 1998h . The broad-band SED of Mrk 421 can often 


be mo deled with a one-zone synchrotron se l f-Compton (SSC ) 
model l lAleksic et al.ll201^ : lAcci^i et alj[201 ll : lAbdo et alj|201 ih 
or with lepto-hadronic models involving proton synchrotron ra- 
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diation and/or photopion interactions (e.g., BQttch er et al.|[2013l : 


iMastichiadis et al.l2013l;lDimitrakoudis et alj 2014). 

Mrk 421 exhibits larg e variations in the TeV, GeV, X-ray, 
and optical wavebands (e.g., Acciari et al.ll201 iLl2009l : iHoran et'H] 
l2009tlBlazeiowski et al.l200 5l). with correlated TeV a nd X-ray vari¬ 
ations (Zg!7|Giebei^^}|20^; l^ssatTet^JI^OO^. The connec¬ 
tion between the TeV and optical bands is less clear and more 
detailed time-dependent mode ls are being deve loped to study the 
co mplex correlations (see e. g. JChen et al.ll2oTll for SSC modeling 
and IMastichiadis et al.ll2013l for lepto-hadronic modeling). 

In 2012 Mrk 421 underwent its largest radio flare ever ob¬ 
served at 15 GHz dHovatta et al.ll2012h . The flare time-scale was 
faster and its amplitude larger than any other radio flare ob¬ 
served from this source in the past 30 years of observations with 
the University of Michigan Radio As tronomy Observatory mon¬ 
itoring program dRichards et al.l 1201^ . The radio flare occurred 
about 40 d after a major flar e was observed in the y-ray band 
dP’Ammando & Orientil 12012) by the Fermi Gamma-Ray Space 
Telescope (hereafter Fermi). Unfortunately Mrk 421 was close to 
the sun at the time of the radio flare so that only limited multi¬ 
wavelength coverage in addition to the radio and y-ray bands ex¬ 
ists, apart from a NuSTAR calibration observation that covered 
the early activity dSalokovic et alj|2013h . In the spring of 2013, 
Mrk 421 underwent anot her major high-energy fl a ring event in the 
X-ray to TeV energies ([Cortina & HoldeJ 1201 3l : iBalokovic et al.l 
I 2 OI 3 I : IPanegue et al]|2013l) . About 60 d later, it was followed by 
fairly small amplitude radio flares in the radio and millimeter bands 
dHovatta et al.ll20l33) . 

In this paper we present the radio data obtained with the 
Owens Valley Radio Observatory (OVRO) 40-m telescope at 
15 GHz, and 95 GHz data obtained with Combined Array for Re¬ 
search in Millimeter-Wave Astronomy (CARMA). We compare the 
radio variations with the y-ray light curves obtained by Fermi, by 
performing a cross-correlation analysis on the full data sets avail¬ 
able from 2008 to 2013. Based on our findings and the coincidence 
of historically rare extreme flares in the two bands, we assume that 
the events in the radio and y-ray bands are physically connected. 
Under this assumption we aim to understand the extreme nature of 
the 2012 radio event by adopting a reasonable physical model with 
the smallest number of free parameters. 

Our paper is organized as follows. We describe the observa¬ 
tions and data reduction in Sect. The light curves and cross¬ 
correlations between the different bands are shown in Sect. [3 We 
model the flares in Sect.|4] discuss our results in Sect.|5] and present 
our conclusions in Sect.|^ Throughout the paper we use a cosmol- 
ogy where Hp = 71 km s“' Mpc^*, Qm = 0.3, and Oa = 0.7 (e.g., 
iKomatsu et al.ll200^ . 


2.1 Radio observations 

Mrk 421 was observed as part o f the blazar monitori ng progranQ 
with the OVRO 40-m telescope dRichards et al.ll20ll]) . In this pro¬ 
gram, a sample of over 1800 AGN are observed twice per week 
at 15 GHz. Mrk 421 has been included since the beginning of the 
monitoring in late 2007. 

The OVRO 40-m telescope uses off-axis dual-beam optics and 
a cryogenic high electron mobility transistor (HEMT) low-noise 
amplifier with a 15.0 GHz center frequency and 3 GHz bandwidth. 
The two sky beams are Dicke switched using the off-source beam 
as a reference, and the source is alternated between the two beams 
in an ON-ON fashion to remove atmospheric and ground contami¬ 
nation. A noise level of approximately 3^ mjy in quadrature with 
about 2 per cent additional uncertainty, mostly due to pointing er¬ 
rors, is achieved in a 70 s integration period. Calibration is achieved 
using a temperature-stable diode noise source to remove receiver 
gain drifts and the flux d ensity scale is deriv ed from observations 
of 3C 286 assuming the iBaars et ^ ( Il977l) value of 3.44 Jy at 
15.0 GHz. The systematic uncertainty of about 5 per cent in the 
flux density scale is not included in the error bars. Complete details 
of the reduction and calibration procedure are given in Richards et 
al. (2011). 

For comparison and to fill gaps in the OVRO sampling, we 
also consider the 37 GHz light curve obtained by the Metsahovi 
Radio Observatory. These observations were made with the 13.7- 
m radome-enclosed telescope using a 1 GHz-bandwidth dual beam 
receiver centered at 36.8 GHz. A detailed descript i on of the data 
reduction and analysis is given in iTerasranta et al.] (Il998h . As the 
uncertainty in the data points is much larger than in the other bands, 
we do not include the data in any subsequent modeling. 


2.2 Millimeter-band observations 

Since February 2013, Mrk 421 was observed 1-3 times per 
week with CARMA as a part of the MARMOtH blazar monitor¬ 
ing project. The observations were made using the eight 3.5 m 
telescopes of the array with a central frequency of 95 GHz and a 
bandwidth of 7.5 GHz. The data were reduced using the MIRIAD 
(Multicha nnel Image Recon struction, Image Analysis and Display) 
software dSault et al.lll995l) . including standard bandpass calibra¬ 
tion on a bright quasar. The amplitude and phase gain calibration 
was done by self-calibrating on Mrk 421. The absolute flux calibra¬ 
tion was determined from a temporally nearby observation (within 
a day) of the planets Mars, Neptune or Uranus, whenever possible. 
Otherwise the quasar 3C 273 was used as a secondary calibrator. 
The estimated absolute calibration uncertainty of 10 per cent is not 
included in the error bars. 


2 OBSERVATIONS 

We present radio and y-ray observations of the two flaring events 
observed in 2012 and 2013. We are especially concerned with the 
extreme 2012 radio event, whereas the 2013 flare is more typical of 
the radio variability observed in Mrk 421. 


2.3 Gamma-ray data 

The y-ray data were obtained with the Large Area Telescope (LAT) 
aboard Fermi, which ob serves the entire sky every 3 hours at ener¬ 
gies of 0.1 - 300 GeV l lAtvrood et al.ll200% . The publicly avail¬ 
able reprocessed Pass 7 datJl were downloaded and analysed us¬ 
ing the Fermi ScienceTools software package version v9r32p5. The 

^ http://www.astro.caltech.edu/ovroblazars/ 

^ Monitoring of y-ray Active galactic nuclei with Radio, Millimeter and 
Optical Telescopes; http://www.astro.caltech.edu/marmot/ 

^ http://fermi. gsfc. nas a. gov/s sc/data/ 
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data w ere binned using the adaptive binning method of iLott et alj 
( l2012h . The uneven bin size was determined in such a way that 
the statistical error in each flux measurement is ~15 per cent 
We used the instrument response functions P7REP_SOURCE_V15, 
Galactic diffuse emission model “glLiem_v05.fits” and isotropic 
background model “iso_source_v05.txt”|3 Source class photons (ev- 
class=2) within 15° of Mrk 421 were selected, with a zenith angle 
cut of 100° and a rocking angle cut of 52°. 

Once the bins were dehned, the photon fluxes in the energy 
range of 0.1 - 200 GeV were calculated using unbinned likelihood 
and the tool gtlike with the Minuit optimizer. All sources within 
15° of Mrk 421 were included in the source model with their spec¬ 
tral parameters, except the flux, fro zen to the values d etermined in 
the 2nd Fermi LAT catalog t2EGL: lNolan et al.]|2012h . Eor sources 
more than 10° from Mrk 421 we also froze the fluxes to the 2EGL 
value. The 10 per cent systematic uncertainty below 100 MeV, 
decreasing linearly in Log(E) to 5 per cent in the range between 
316 MeV and 10 GeV and increasing linea rly in Log(E) up to 15 
per cent at 1 TeV jAckermann et alj|2012h is not included in the 
error bars. For the purpose of the modeling presented in Sect. [4] we 
also obtained the energy fluxes, with a power-law spectral model 
where the index is frozen to the 2FGL catalog value of T = 1.77. 


3 LIGHT CURVES 

The radio, millimeter and y-ray light curves are shown in Fig. [T] 
They cover the time range since the beginning of the Fermi mis¬ 
sion in August 2008 (MJD 54688) until the end of October 2013 
(MJD 56610). The light curves illustrate the unusual nature of the 
2012 flare (flare 1) in both radio and y rays, that lasted from May 
2012 until October 2012 (MJD 56060 - 56225). This is a unique 
event, especially in the radio band where such fast and promi¬ 
nent flares have not been observed before in Mrk 421. The 2013 
flare (flare 2) from March 2013 until November 2013 (MJD 56350- 
56610), although prominent, is much broader and of lower ampli¬ 
tude. 

A full cross-correlation analysi s between more than four year s 
of OVRO and LAT data was done bv iMax-Moerbeck et alj l l2014bh . 
They used weekly binned y-ray light curves and data from 2008 
until November 2012, including the rapid 2012 flare. They found a 
peak in the discrete correlation function (DCF), with the y rays 
leading the radio by 40 ± 9 d. The significance of this corre¬ 
lation was between 96.16 and 99.99 per cent depending on the 
power spectral density (PSD) model used for the light curves, with 
a best-fitting value o f 98.96 per cent. For details of the signifi¬ 
cance estimation, see iMax-Moerbeck et alj l l2014ah . We repeated 
the cross-correlation analysis using the extended light curves con¬ 
sidered here, and found that the DCF shows a broad peak (~ 30 d). 
In particular, the time delay ranges from 4 0 to 70 d with y rays 
leading , consistent with the estimate from iMax-Moerbeck et al.l 
( l2014hl) . The significance of the peak was from 91.90 to 99.99 
per cent, with a best fit value of 97.36 per cent, depending on the 
PSD model. The difference co mpared to the exact value derived by 
iMax-Moerbeck et al.l ( l2014bl) is because they only considered the 
maximum of the DCF, while the peak itse lf is bro ad. 

Recently, lEmmanoulopoulos et al.l l l2013h introduced a 
method for simulating light curves, which also accounts for their 
flux distribution. This is more appropriate in the y-ray light curves. 


^ http://fermi.gsfc.nasa.gov/ssc/data/access/lat/BackgroundModels.html 
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Figure 1. From top to bottom: light curves of Mrk 421 in y rays from Fermi 
LAT, radio 15 GHz from OVRO, 37 GHz from Metsahovi, and 95 GHz from 
CARMA. The light curves span a time period from August 2008 until 
November 2013, except for the CARMA light curve that starts in Febru¬ 
ary 2013. The solid lines indicate the time range used to model the 2012 
flare from 13 May 2012 (MJD 56060) until 25 October 2012 (MJD 56225). 
The dashed lines indicate the time range for the 2013 flare from 9 March 
2013 (MJD 56360) until 14 November 2013 (MJD 56610). Because of the 
larger uncertainties in the data, the 37 GHz light curve is not used in any 
subsequent modeling. 


as the light curves have a non-Gaussian photon flux distribution. 
Using this method, the significance of the correlation increases to 
99.82 per cent, when using the best-fitting PSD. In the rest of the 
paper, we will thus assume that the major flares in radio and y rays 
are physically connected. 

The y-ray flares in both 2012 and 2013 have significant sub¬ 
structure, as already shown in Fig.[T](top panel). If we also consider 
the broadness of the DCF peak, it becomes unclear which y-ray 
spike of the overall flare is most plausibly associated with the radio 
peaks. Therefore, we consider both alternatives in our modeling. 


3.1 The flare of 2012 

On 16 July 2012 (MJD 56124) the y-ray flux reached the highest 
value since the start of the Fermi mission jD’Ammando & Orientil 
l2012l) . The y-ray flux increased by a factor of three from (3.4 ± 
0 .6)x 10“’ phcm“’ s“* to (Ll±0.2)x 10“® phcm“’ s“* within seven 
days. The flare appears double peaked with the second flare peak¬ 
ing on 15 August 2012 (MJD 56154) at a flux of (9.5 ± 0.1) x 
10“’ phcm“’ s“'. The total duration of the flaring event was about 
60 d. The OVRO 15 GHz light curve exhibits a fast rise which leads 
to an increase of the flux density by a factor of about two, i.e. 
from 0.6 to 1.1 Jy, in the period 6 August 2012 - 21 September 
2012. This flux density is higher than any flux density measured 
at 15 GHz in the OVRO program or during the preceding 30-plus 
years of monitorin g with the University of Michigan Radio Astron¬ 
omy Observatory i Richards et alJlTofsh . 

At the time of the first y-ray peak there is a gap in the OVRO 
15 GHz light curve. Thus, a double peaked radio flare cannot be 
excluded a priori. However, a higher frequency radio light curve 
observed at 37 GHz at Metsahovi Radio Observatory (third panel 
from the top in Fig.|T]( shows only a single flare during 2012. Given 
the spectral proximity of the two radio bands and the similar fea- 
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Figure 2. y-ray flux (top) and 15 GHz flux density (middle) during the 
2012 flare. The discrete con'elation function between the two light curves 
is shown in the bottom panel. A negative time delay means that y rays are 
leading the radio. The (blue) solid line shows the exponential fit that is used 
to estimate the variability brightness temperature and Doppler factor. The 
residuals of the fit are on average less than 20 mJy. The time range in the 
figure corresponds to the solid lines in Fig.[T] 


tures in the two light curves, it is safe to assume that the OVRO 
sampling does not miss a peak at the beginning of the event. Be¬ 
cause of the larger statistical uncertainties in the 37 GHz data, we 
do not include them in our subsequent analysis as they would not 
add additional constraints on the modeling. 

In order to obtain a better estimate of the time delay between 
the y-ray and 15 GHz light curves for this flare only, we use the 
DCF method over the time period of the flare ( Fig. [T] bottom 
panel). There are two peaks in the DCF which simply correspond 
to the delays between each of the spikes in the double-peaked y-ray 
flare and the single radio flare. The DCF peaks are both consistent 
with the delay measurement we obtained for the full light curves. 
Although the correlation is stronger for the ~40 d lag, the ampli¬ 
tude of the DCF peak at about -70 d is not low enough to justify 
exclusion of a possible association of the radio flare with first y-ray 
spike. Thus, we will test both possibilities in Sect.|4] 

We can estimate the Doppler boosting factor of the radio flare 
by assuming that the rise time of the flar e corresponds to the light 
travel time across the em ission region jLahteenmaki & Valtaoial 
1 19991 : iHovatta et alj|2009l) . We fit an exponential function of the 
form 

S(0 = AS (1) 

to the light curve, where AS is the amplitude of the flare, t,„ax 
is the peak location of the flare and Cue is the rise time of the 
flare. The dec ay time of the flare has been fro zen to 1.3 times the 
rise time, as in lLahteenmaki & Valtaoial l ll999l) . The fitting is done 
using the MultiNest Markov-Chain Monte Car l o (MC MC) algo¬ 
rithm jFeroz & HobsonlEoOSl i lFeroz et al.|[200^l2013h . As shown 
in Fig. 1^ middle panel, an exponential function fits the 15 GHz 


flare fairly well. From the fit we obtain the rise time of the flare 
trise = 10-6 ± 0.5 d and amplitude AS = 0.52 ± 0.01 Jy, where the 
uncertainties are estimated from the MCMC analysis. 

Assuming an emission region with the geometry of a uniform 
disk, we can estimate the lower limit of the variability br ightness 
temperature jLahteenmaki et al.lll999l : lHovatta et alj2009l) 


Tva.- = 1.548 X 10^ 


ASrff 


+z)’ 


( 2 ) 


where is the luminosity distance to the object in meters (here, 
4.0 X 10^"* m), z is the redshift, v is the frequency in GHz, t,ise is in 
days, and AS is in janskys. This gives us - (5.2 ±0.5) x 10'^K. 

The variability brightness temperature is related to the 
Doppler factor, 6,^s 6 = (rv ar/rb.int)‘^^, where TKim i s the intrinsic 
brightness temperature (e.g., Lahteenmaki et alJl999h . As the vari¬ 
ability brightness temperature estimate is a lower limit, the Doppler 
factor estimates are also lower limits. The largest uncertainty in the 
Doppler factor estimate comes from the uncertainty in the intrin¬ 
sic brightness temperature and what method is used to estimate it. 
Assuming equi partition between the particles and magnetic field, 
Tb.int = 10' ‘ K jReadhead 1 1 994h . This results in a Doppler factor 
6 > 3J. If we use a value of Tbint = 5 x 10*“ K, as determined by 
iLahteenmaki et all ( Il999h . the Doppler factor is <5 > 4.7. 

We can further constrain the intrinsic brightness tempera¬ 
ture by comparing the variability brightness temperature with 
the brightness temperature obtained via simultaneous Very Lon. 


_ TV J _ 

Baseline Array (VL BA) observations jLahteenmaki et aH 1 19991 : 
iHovatta et alj2013^ . In order to study the parsec-scale jet structure 
after the 2012 flare we conducted 5 epochs of target-of-opportunity 
VLB A observations of Mrk 421 at several frequency bands (J. 
Richards et al. in preparation). The first one of these was taken on 
12 October 2012 when the radio flare was already decaying. The 
brightness temperature estimate from the 15 GHz data is Tvlbi = 
5.2 X 10*“ K (for a uniform disk), which depends on the Doppler 
factor as (5 = TvLBi/Fbjnt- We can then solve for the intrinsic bright¬ 


ness temperature by calculating Tbjnt = “ 5.2x 10“ K. 

If we use this value for the intrinsic brightness temperature and the 
variability brightness temperature, we obtain 6 > 10. Thus, we con¬ 
clude that the lower limit of the Doppler factor is <5 ~ 3 - 10. 

We note that the intrinsic brightness temperature obtained us¬ 
ing the VLB A data is about 10-20 times below the equiparti¬ 
tion limit. We think the peak brightness temperature is likely higher 
than our estimate because by the time of the first VLBA epoch, the 
single-dish flux density had already declined by 30 per cent from 
the peak. Therefore it is likely that the true simultaneous bright¬ 
ness temperature from the VLBA observations is at least 30 per 
cent higher, because the core could have also been more compact, 
increasing the brightness temperature even further. Because of the 
strong dependence of Tb.int on Tvlbi, any uncertainties in the VLBA 
parameters are magnified in the estimate of Tb.mt- A slightly h igher 
Tvlbi would also agree with estimates from iLico et al.l j2012l) . who 
found the core brightness temperature of Mrk 421 to be of the order 
of few times 10** K, in agreement with equipartition arguments. 


3.2 The flare of 2013 


In Apr il 2013, Mrk 421 was again flaring i n the X-ray to TeV 
bands jBalokovic et al.l 1201 3l : IPanegue et al.l 12013 ), reaching the 
highe st levels ever observed at TeV energies Jcortina & HoldeJ 
l2013h . Triggered by this activity, we began monitoring the source 
more frequently at CARMA. 
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The appearance of this flare was very different from the 2012 
flare, both in the y rays and radio. In y rays the activity began in 
March 2013 reaching the highest peak on 14 April 2013 (MJD 
56397). The flux increased from about 2 x 10“’phcm^^ s“* to 
(6.7 ± 0.9) X 10“’ ph cm“’ s“* over about 30 d and the flaring period 
consisted of several peaks over a total duration of about 100 d. At 
15 GHz, the flux density began increasing in April 2013 from about 
0.6 Jy to its peak of 0.77 Jy on 18 June 2013 (MJD 56461). Sim¬ 
ilarly, at 95 GHz the flux density increased from about 0.5 Jy to a 
peak of 0.71 Jy on 8 June 2013 (MJD 56451). 

The DCFs between all the bands are shown in Fig.[3 The high¬ 
est correlation is found for the 15 and 95 GHz data, with a peak at a 
delay of about -10 to -20 d, with 95 GHz leading. The time delays 
between the CARMA (OVRO) and LAT data can only be estimated 
at about -60 d (-70 d) because of the broad peaks of the DCF. 
These estimates are, however, compatible with the broad peak ob¬ 
tained from our cross-correlation analysis of the full light curves, 
and thus we assume that the events are physically connected in the 
radio and y-ray bands. 


4 CONDITIONS IN THE FLARING REGION 

In this section we attempt to explain the rough features of the 2012 
and 2013 flaring periods, e.g., time delays between various energy 
bands, peak fluxes and pulse profiles, by adopting the simplest pos¬ 
sible theoretical framework (the single zone SSC model), together 
with a minimum set of free parameters. In particular, we restrict our 
modeling to the major flares observed in the periods of MJD 56060- 
56225 and MJD 56360-56610, which are denoted as “flare 1” and 
“flare 2”' in Sect.[3(see Fig.[TJ, while modeling of the smaller am¬ 
plitude variability seen in both bands lies out of the scope of this 
work. As the ‘goodness’ of the fits was not the main point of inter¬ 
est here, we did not attempt a detailed parameter space search for 
finding the set with the best^’ value. 

Based on the long-term radio light curve shown in Fig.[T] the 
2013 radio and millimeter-band flares are more typical of blazar 
radio emission than the 2012 radio flare: they are less sharp, have 
a decay timescale larger than their rise timescale, and the flux den¬ 
sity increases by no more than a factor of ~ 1.4 in both the radio 
and millimeter bands. We start our analysis with the 2013 flaring 
events in the context of a typical SSC model. Then, we highlight 
the differences between the 2012 and 2013 radio flares and discuss 
possible modifications within the same framework that may explain 
the 2012 data. Our main goal is to investigate what conditions are 
required to produce this extreme and unique radio flare. As already 
discussed in Sect.[2 we assume that the events in the radio and y- 
ray bands are physically connected. We note that, in what follows, 
we ignore the small differences between quantities measured in the 
observer frame and the source frame since 1 -l- z 1. 

4.1 The y-ray and radio flares of 2013 

Motivated by the fact that the CARMA and OVRO fluxes peak ~ 60 
and 70 d after the y-ray flare, respectively (Sect. 13.2b . and by the 
fact that the decaying part of their light curve is approximately ex¬ 
ponential, we applied the simple scenario where both the y-ray and 
radio flares originate in the same region, which does not necessarily 
have the same physical conditions as the region responsible for the 
quiescent emission. We note that we do not attempt to model the 
sub-structure of the flares, but we only consider the highest peak in 
each band. In particular, we focus on the largest y-ray flare, which 





Time (MJD) 

Figure 4. Fermi-LAT light curve (top), CARMA 95 GHz light curve (mid¬ 
dle) and OVRO 15 GHz light curve (bottom). The results of the one-zone 
SSC model described in the text are shown with blue lines in ah panels. The 
time range of the ht corresponds to the 2013 flare indicated by the dashed 
lines in Fig.[T] 


peaks at MJD ~56400. We also note that the DCF only shows a 
single peak for this flare, unlike in 2012 where we will consider 
multiple associations between flares. 

In our scenario the y-ray flare is produced by an instantaneous 
injection event of electrons having Lorentz factor yo and emitting in 
y rays through Compton scattering of synchrotron photons (SSC). 
In particular, electrons are injected at t = 0, which is set to be 
the time of the y-ray flare, and then they are left to evolve via 
synchrotron and SSC cooling. In this context, the observed delays 
fh ~ 60 d and ff ~ 70 d between the y-ray and the radio flares at 
Vh = 95 GHz and Vf = 15 GHz, respectively, correspond to the 
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Figure 3. Top: Light curves at y rays (left), 95 GHz (middle), and 15 GHz (right) during the 2013 flare. Bottom: Discrete con'elation functions between the 
light curves. Negative time delay means that higher frequency emission (y-ray or 95 GHz) is leading. The time range corresponds to the dashed lines in Fig.[T] 


cooling timescale of electrons that have been injected at t = 0 with 
Lorentz factor yo- 


4.1.1 Analytical estimates 


To define the size R of the emitting region we choose as a typ¬ 
ical variability timescale (fy) the one dictated by the y-ray light 
curve. For the purpose of our analysis, we choose a variability 
timescale of ty = 20 d based on the light curve. We note that 
the exact value is not critical as the estimate will be refined in the 
next section where we model the flares numerically. We find that 
R x 2x 10 *’ cm (6/4) {t.,/20 d), where we normalized the Doppler 
factor to 4 (see Sect. [IB- In this framework, one can derive the 
magnetic field strength B of the emission region as a function of 
the radio frequency Vh, the observed time delay 4 and the Doppler 
factor 6. Electrons that have been injected at t = 0 with Lorentz 
factor 7 o cool due to synchrotron losses and reach a Lorentz fac¬ 
tor 71 ,. This can be found by solving the characteristic equation for 
synchrotron cooling and is given by 


7h = 


1 

1/70 + a*h 


(adth) * , 


(3) 


where a = (4/3)(TTc(7B/meC^ and Ub = B^/&n is the energy den¬ 
sity of the magnetic field. The approximation holds as long as 
Jo 7 h. Combining Eq. © with the characteristic synchrotron 
frequency vj, = 6m^c^bj^/h, where b = B/B^, 6^ = 4.4 x 10*’ G, 
and h = 6.63 x 10“” erg s is the Planck constant, we find that the 
required magnetic field strength of the region is 


B^O.SG 



which depends most strongly on the time delay between the 7 -ray 
and radio flares. Notice also that lower values of the Doppler factor 
favour stronger magnetic fields. The expected time delays for the 
two frequencies should satisfy the relation 4/4 = (vh/vf)*’’. This 
is roughly consistent with the observed values, since VITty = 2.5 
and 4/4 is in the range of 1 - 1.7 based on the broad peaks in the 
time delays (see Fig.©. The duration of a flare is also related to the 
magnetic field strength, since it depends on the electron cooling 
timescale. For a given frequency Vh oc By’ the full width at half 
maximum (FWHM) of a flare will scale roughly as oc oc 

Thus, for strong enough magnetic fields the duration of 
a flare even at low frequencies may be shortened significantly. 

Assuming that the inverse Compton scattering occurs in 
the Thomson limit, which will be checked a posteriori, the 
typical energy of up-scattered photons is given by Ey 
(4/3)6m,.c^byg. Using the value of B derived above, we find 
that the injection Lorentz factor of electrons is 70 =: 2 x 

10'* (Ey/4 GeV)‘^'' (6/4)-''* (B/0.5 G)-''*, where we chose Ey = 
4 GeV as a representative value for the energy of y-ray photons. 
Because of the weak dependence of 70 on Ey, a different choice 
of the y-ray energy would not affect the derived value for the in¬ 
jection Lorentz factor. Increasing Ey by a factor of 50, for ex¬ 
ample, would increase 70 by only a factor of 2. In principle, 
electrons could have been injected with 7 > 70 and still produce 
the GeV flare, since they would very quickly cool to the value 
70 . In this case a contemporaneous TeV flare would also be ex¬ 
pected, and indeed one was seen by the MAGIC and VERITAS 
instruments dCortina & Holdeiil20L3h . Having derived the expres¬ 
sion for 7 o, which does not depend strongly on the parameters, we 
can now verify our earlier assumption that the Thomson limit ap- 
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plies. Since byl = 0.09 (f) < 314, we are 

safely in the Thomson regime. The derived values for the mag¬ 
netic field strength and the Lorentz factor of injected electrons are 
reasonable and consistent with typical SSC models of high fre¬ 
quency peaked BL Lac ob jects, such as Mrk421 fe.g.. lAbdo et al.l 
I 2 OI ll i fAleksic et ^l2012h . Note, however, that here we adopted a 
low value for the Doppler factor implied by radio observations, in 
contrast to typical SED modeling where much larger v alues, e.g. 
20 - 50, are usually used (see e.g.. lMaraschi et aljl999ll 

4.1.2 Numerical results 

The required electron injection luminosity is determined numer¬ 
ically by fitting the observed energ y fluxes. For this, we em¬ 
ploye d the numerical code described in iMastichiadis & Kirkl ( ll91^ 
Il997h . Using as a stepping stone the values determined analyti¬ 
cally in Sect. 4.1.1., we derive the following set of parameters: 
R = lO'^ cm, 6 = 0.1 G, (5 = 2.2, 70 = 2 x 10“^ and the elec¬ 
tron injection compactness = 1.6 x 10“^. This is defined as 
= (TTRUelm^c- where 11^ is the energy density of electrons at 
injection time as measured in the rest frame of the emission re¬ 
gion and CTj- = 6.65 x 10“^^ cm^ is the Thomson cross section. 
Although our analytical estimates were based upon the hypothesis 
of instantaneous injection, here we find that an injection episode 
lasting ~ Rjc, or equivalently 20 d, is better in reproducing the 
observed flares. 

Our results are illustrated in Fig.|4l where we use energy fluxes 
to facilitate comparison between the diiferent bands. In the radio 
band, flux densities are converted to energy fluxes by multiplying 
by the observing frequency and converting to cgs units. The en¬ 
ergy fluxes in the y-ray band are obtained as described in Sect. 2. 
The model light curves can describe the radio observations fairly 
well, although the model y-ray flare is slightly broader than the ac¬ 
tual data. Still, the results of Fig.|4]are satisfactory given the small 
number of free parameters and the fact that we have not attempted 
to find the best fit (i.e., with the lowestvalue). The radio flares in 
both frequencies have wide pulse profiles and exhibit a slow expo¬ 
nential decay after their peak. Both of these features contrast with 
the sharpness and the symmetry of the 2012 flare. In the present 
scenario, the width of each flare is related to the cooling timescale 
of electrons emitting at the particular energy band and, thus, the 
radio flares are wider than the one in y rays (see Sect. 4.1.1.). The 
asymmetry of the pulse profiles at lower energies is a strong predic¬ 
tion of this model. Possible expansion of the flaring region and/or 
decay of the magnetic field would increase the predicted asymme¬ 
try. 

For the derived parameters, the region is initialljH particle 
dominated with t/e ~ 0.18 erg cm“^ and Ub = 4 x erg cm“^. 
The value of the Doppler factor derived here is approximately half 
of that obtained in Sect. [3T] assuming equipartition. However, the 
large ratio UsIUb - 450 derived from the values menti oned above 
would also imply a higher brightness temperature (e.g.. lReadheadl 
Il994h by a factor of 450*^* ~ 2 and thus a lower value of S by 
a factor of 2“'^^ ~ 0.7. Given that the fit shown in Fig.|4]is not 
unique, one could find parameters that would decrease the initial 
ratio UcIUb and bring the emission region closer to an equiparti¬ 
tion state. By assuming that a fraction 77rad of the injected luminosity 
in electrons is radiated as y rays, i.e. t/rad^e"'' = Lf’‘crjl4nRS*m^c^, 

^ The electron energy density will actually decrease with time because of 
(i) cooling and (ii) no replenishment of particles in the region. 


and using the definition of as well as the equations for R and B 
in Sect. 4.1.1., we find that 



where Fq = Ffyi x 10“'* erg cm“^ s“*. The ratio of energy den¬ 
sities depends strongly on the Doppler factor. Thus, searching for 
possible fits with a slightly higher value of the Doppler factor, e.g., 
6 — 1, would bring down the ratio from 450 to close to unity and 
would ensure rough equipartition between the particles and mag¬ 
netic field at injection. 

4.2 The y-ray and radio flares of 2012 

The 2012 radio flare is a unique event not only because of its large 
flux increase but also because of its symmetric pulse profile, which 
resembles flares at higher energies where the corresponding radi¬ 
ating particles cool efficiently. According to the cross-correlation 
analysis presented in Sect. 3.1, the 2012 radio flare may be asso¬ 
ciated with either the first or second spike of the y-ray flare (see 
Figs. 1 and 2). We investigate both possibilities under the assump¬ 
tion that there is indeed a physical connection between the extreme 
radio and y-ray flares. For this we begin with the parameters we de¬ 
rived for the 2013 “typical” flare, and examine two alternative sce¬ 
narios. These include changes of the Doppler factor and magnetic 
field strength, since both parameters are related to the time delay 
(see Eq. II and to the pulse profile (see Sect. 4.1.1.). We emphasize 
that in both scenarios we attempt to introduce minimal changes to 
the parameters. As before, we do not consider the small-amplitude 
sub-structure of the flares. 


4.2.1 Association of the radio flare with the first y-ray spike 

We first investigate the possibility that the radio flare is associ¬ 
ated with the narrower first spike of the y-ray flare, which peaks 
at ~56133 MID. In order to do this, we use the same parameter 
set as derived for the 2013 flare, but add variations in the Doppler 
beaming to the model. Such changes could occur if the emission 
region moves on a curved trajectory with a changing viewing an¬ 
gle, or if the bulk Lorentz factor of the emission region changes. 
Varying Doppler factors have been previously us ed to explain radi o 
spectral variations in, for example, S5 0716+714 ( iRani et alj2013h . 
The Doppler factor 6 is modeled as 

6{t) = do (1 + g(r) cofllmlP)), ( 6 ) 

where do = 2 . 1 , r is the time in the comoving frame in R/c unitsQ, 
P = 5.4 is the period of the variation in Rjc units and g(T) is the 
function 

g(T) = 0.3, for T < Tbr (7) 

g(T) = forr Tbr, ( 8 ) 

with Tbr = 1-2R and T = Q.2P. 

The model-derived pulse profiles are shown in Fig.|5]with blue 
solid lines. Although the cooling timescale of electrons emitting at 
radio frequencies is long compared to those emitting initially at y 
rays, the shape of the radio flare (right panel in Fig.O is well re¬ 
produced here, because of the change in the Doppler factor, which 

® Time is measured with respect to the time of the first spike seen in the 
2012 y-ray flare. 
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Figure 5. Fermi-LAT light curve (left) and OVRO 15 GHz light curve (right) in energy flux units along with the light curves (blue lines) obtained from the 
model. Solid and dashed lines demonstrate the model results, if the radio flare is associated with the first and second spikes of the y-ray flare, respectively. A 
different x-axis range is used in the right panel to show the extended model light curve more clearly. 


does not allow the observer to see the typical wide pulse profile (see 
e.g., Fig.H. The effect of the variable 6 on the y-ray light curve is 
negligible because the radiating electrons at the time of injection 
have short cooling timescales. Note that the exponentially decay¬ 
ing amplitude in the variation of 6 is required in order to avoid any 
excess of the radio emission following the extreme flare. 

Assuming that the change in the Doppler factor is caused by 
small variations in the viewing angle 0 while the Lorentz factor 
remains constant at T = 10, no extreme variations in either 6 or the 
angle 9 = cos^' (l - are required^In particular, 6 and 

9 change by less than 28 and 22 per cent, respectively, relative to 
their average values. 

The fact that the necessary parameter changes are small make 
this scenario plausible and attractive. However, if we attempt to in¬ 
terpret the second spike of the y-ray flare in the same framework, 
then we cannot avoid the appearance of a radio flare ~ 70 d later, 
i.e. at MJD 56220. Its absence implies one (or more) of the follow¬ 
ing: (i) a faster onset of the exponential decay of d; (ii) a different 
functional form for 6{t)\ (iii) different conditions in the emission 
region which suppress the 15 GHz emission, such as different mag¬ 
netic field and/or size; (iv) a different framework for the second 
y-ray spike. 

Finally, the choice of a damped oscillator for modeling 6 {t), 
albeit plausible, was not based on a physical picture. According to 
the discussion of Sect. 4.1.1., the FWHM of a flare depends on the 
magnetic field strength. Therefore, another alternative might be to 
adopt a stronger magnetic field for the emission region. Although 
a stronger magnetic field would lead to a shorter duration of the 
flare, an even lower value of the Doppler factor would be required 
in order to explain the observed time delay between the y-ray and 
radio flares (see Eq.|4l(. Given that the value of the Doppler factor 
already lies at the low end of values consistent with observations, 
we do not examine this possibility in more detail. 


^ This Lorentz factor is approximately five times higher than the value de¬ 
rived in Sect.[^for the bulk Lorentz factor of the jet (Lyar = 1.9). If we had 
used this value instead, only the angle 9 would be affected, i.e. it would be 
larger by a factor of 2. However, it is still possible that the Lorentz factor of 
the em ission region is higher than that of the jet (see e.g., lNalewaiko et*^ 

l20il . 


4.2.2 Association of the radio flare with the second y-ray spike 

We now consider the scenario where the radio flare is associated 
with the second spike of the y-ray flare. This scenario would arise 
naturally if the first y-ray spike occurred when the emission re¬ 
gion was still optically thick to radio emission. In this case, the 
time-delay between the y-ray and radio flares is shorter than be¬ 
fore, which, in our framework, suggests faster electron cooling. We 
searched, therefore, for reasonable fits using a stronger magnetic 
field than that adopted in Sect. 4.2.1, and obtained the following 
parameter set: B = 0.25 G, d = 2.3 and = 1.3 x 10“^. All other 
parameters, including the injection profile, are the same as before. 

The resulting pulse profiles are shown in Fig. with dashed 
blue lines. We find that the light curves obtained by the model 
are in rough agreement with the observations, with the radio light 
curve having a slightly longer decay timescale than is observed. In 
this scenario, even without Doppler factor variations, we obtain a 
sharper radio flare than the one in 2013 (see e.g. Fig. 4). Compar¬ 
ing these parameter values with the first scenario, we see that al¬ 
though we adjusted the numerical values for three parameters, only 
the magnetic field is significantly altered. It has increased by a fac¬ 
tor of 2.5. We can compare the FWHM of the 15 GHz radio flares 
in 2012 and 2013 as obtained by our model (Figs. |4] and |5j. We 
find their ratio to be AT 2012 !A.T 2013 . ~ 0.16, where AT denotes the 
duration at FWHM. This is comparable to the analytical estimate 
we can obtain from Eq. El A7’20l2/Ar2013 - (S 2 OI 2 /S 2 OI 3 ) - 

(0.25/0.1)“^/^ = 0.25. 

We note that this model is subject to tight constraints on its 
parameters. These must be adjusted to ensure: (i) the electron cool¬ 
ing timescale results in the observed delay between the y-ray spike 
and the radio flare (see Eq. 0); (ii) the cooling timescale of elec¬ 
trons radiating at the radio frequency is such that it explains the ob¬ 
served width of the pulse; and (iii) the synchrotron self-absorption 
frequency, which depends on , B and R, is below the radio fre¬ 
quency where the extreme flare is observed. From the above it be¬ 
comes clear that fine tuning of the parameters is required for this 
scenario to be viable. 


5 DISCUSSION 

Unlike in many other blazars, prominent radio flares are rare in 
Mrk 421 and there are few examples where the radio flares have 
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been modeled in detail. In 1997 a 22 GHz flare wa s observed 
with the Metsahovi Radio Observatory 14-m telescope jTosti et al.l 
1 19981) . The flare lasted about 200 d during which the flux den¬ 
sity increased by a factor of two. The radio flare occurred 60 d 
after a large optical flare that also included very high optical po¬ 
larization. Due to the high optical polarization during the flare, 
the authors interpreted t he flare within the shock-in-jet model of 
iMarscher & Gead ( Il985h . The delay between the optical and radio 
bands was explained by opacity effects. While this model may ex¬ 
plain the general features of the flares and the delay between opti- 
cal and radio ban ds, more detailed modeling was not presented by 
iTosti et al] ( ll998h . 

In 2001 a sim ultaneous flare was obse rved at TeV, X-ray and 
5 GHz radio bands iKatarzvhski et al.l2003h . The radio flux density 
increased by a factor of 1.5 and the flare lasted only for 15 d. The 
absence of time delays between the bands indicates that the emis¬ 
sion region was optically thin even at radio frequencies. An instant 
injection of particles into the jet base was unable to explain the ra¬ 
dio flares as the electrons would have cooled too much to generate 
significant flux density changes before reaching the optically thin 
regime for radio emission. Instead, detailed modeling of the flare 
suggested that it was caused by in-situ acceleration of electrons in 
a blob in the jet. The expansion of the blob explains the decay of 
the flare at all bands. 

The 2013 flaring event was a more typical example of radio 
flares in Mrk 421, with a broad flare in the millimeter and cen¬ 
timeter bands (see Fig. [T] for the long-term behavior). Assuming 
that the events in the different bands are physically connected, we 
modeled the most prominent features of this event with a simple 
one-zone SSC model where an instantaneous injection of elec¬ 
trons produces a y-ray flare. Then, as electrons cool because of 
synchrotron and SSC energy losses, they radiate at longer wave¬ 
lengths and eventually produce the radio flares with a delay rel¬ 
ative to the y-ray flare. Within this context, we estimated the re¬ 
quired magnetic field strength, the size of the emission region, and 
the Lorentz factor of electrons at injection. While the simple model 
presented here can describe the main features of the data satisfac¬ 
torily, more information from simultaneous multi-wavelength ob¬ 
servations is needed in order to discriminate between other pos¬ 
sible models. For example, one could postulate a scenario where 
the delay between the y-ray and radio flares is related to opac¬ 
ity effects, such as if the source were optically thick in radio fre¬ 
quencies at the time of the GeV flare but eventually became opti¬ 
cally thin due to the exp ansion of the emission region or decay of 
the magnetic field (e.g.. lFuhrmann et alj|2014l) . Although the co¬ 
incident unprecedented events between the LAT and OVRO data 
suggests a physical connection of the high- and low-energy emis¬ 
sion, models where two spatially separated emission regions are in¬ 
voked l lBlazeiowski et al][2005l : IPetropouloull2014h are still viable 
alternatives. We note that the 2013 flare occurred during a planned 
multi-wavelength campaign and Mrk 421 was observed regularly 
with numerous instruments at various frequency bands. Therefore 
we anticipate several dedicated studies of the behavior of Mrk 421 
during this flare. 

Even though the 2012 radio flare was more extreme, the 
Doppler factor inferred from the radio variability time-scale is 
fairly low, only about 3-10. This is in accordance with the low ob¬ 
served apparent speeds obtained through radio interferometric ob- 
servations where only subluminal speeds have been detected (e.g., 
IPiner et alJl999l : [Lico et af]|2012l : lLister et alj|2013h . Based on the 
jet brightness asymmetry, low observed appare nt speeds, and low 
brightness temperature of the core component, iLico et ti] |20l^ 


estimate the radio Doppler beaming factor to be about 3, consistent 
with our estimate. If simultaneous X-ray or TeV data were avail¬ 
able, we could estimate a minimum Doppler factor by demand¬ 
ing the sour ce to be optically thin to the highest energy emitted 
photon (e.g.. lDondi & Ghiselliniiri995h . From previous TeV flares 
of Mrk 421 it has been estimated that the minimum Doppler fac- 
tor required for th e TeV emission to be optically thin is ~ 10 
dCelottietalJI 19981) . which is marginally higher than the value ob¬ 
tained from our radio data when equipartition is assumed. 

Preliminary analysis of our follow-up VLBA observations in¬ 
dicate that the component speeds were consistent with subluminal 
motion even after the major flare dRichards et al.ll20l^ . If we take 
the fastest observed jet speed, /t^pp = 0.28c, ob tained from severa l 
years of monitoring at 15 GHz with the VLBA dLister et alj|2013h . 
we can estimate the jet Lorentz factor as 




(9) 


where dvar is the Doppler factor inferred from the variability time- 
scales. This results in a Lorentz factor of Fvar ~ L7 - 5. Similarly, 
we can estimate the viewing angle of the jet to be 


Avar = arctan ■ 


2j8ap 


( 10 ) 


resulting in 6var ~ 0.3 - 4°. As noted by iLico et al.l d2012h . it is 
unlikely that the component speeds resemble the flow speed in 
Mrk 421 because unreasonably small viewing angles would be re¬ 
quired to explain the beaming estimates from the jet/counter-jet ra¬ 
tio. In this case the above equations would not be valid for estimat¬ 
ing the true Lorentz factor and viewing angles. However, they find 
a viable scenario for Mrk 421 with a structured jet, where the view¬ 
ing angle is between 2° and 5°, and the Lorentz factor of the radio 
emitting region about 1.8, in agreement with our estimates. 

One of the long-standing problems in modeling of the high 
synchrotron peaked blazars is the large discrepancy between the 
Doppler factor values inferred from radio observations and those 
obtai ned by SEP mod e ling, with the form er being usually < 10 
(e.g.. IPiner et al.l [79991 : iLister et alj 2013) and the latter lying i n 
the range ~ 20 - 50 (e.g.. lMaraschi et ail 19991 : lAbdo et al.ll201 ll) . 
Several alternatives have been explored to explain this well-known 
“Dop pler factor crisis”, s uch as a structured jet i Ghisellini et al.l 
1200^ . a decelera ting jet dG eorgano poulos & Kazanasll2003h . aiid 
a jet-in-jet model jGiannios et alj 20091) . In the structured jet model 
the high-energy emission comes from a fast spine of the jet while 
the ra dio emission is produced in a slower sheath jGhisellini et al.l 
l2005h . This model is favored by observations of limb brightening 
of the Mrk 421 jet at 43 GHz JPiner et al.ll2010h . 

In the present work, where we have adopted a single-zone 
emission model, we attempted to avoid a large discrepancy between 
the Doppler factor values inferred from the observations and those 
used in the modeling of the flares. However, one could relax this 
condition and search for possible fits to both the radio and y-ray 
flares with S > 20. One can estimate the effect of a higher 6 on 
our results by inspection of the analytical expressions derived in 
Sect. 4.1.1. In particular, we find that R oc S, B oc jq oc 
UcIUb and FWHMoc 6'^^, where we assumed that the ra¬ 

dio observing frequencies and their time-delays with respect to the 
y-ray flare are fixed. On the one hand, choice of a higher d would re¬ 
quire a weaker magnetic field, a larger emission region and would 
reduce the ratio of particle to magnetic energy densities (see dis¬ 
cussion in Sect. 4.1.2). On the other hand, the model light curves 
would be wider than those presented in Fig. 4, since the FWHM 


© 2014 RAS, MNRAS 000. [11011 




























































10 T. Hovatta et al. 


would increase. This is a direct result of the longer electron cooling 
timescale due to the weaker magnetic field. We conclude, there¬ 
fore, that values > 20 are less plausible for the modeling of both 
the radio and y-ray flares. Another possibility is to assume a high 
Doppler factor value for the y-ray flare alone and a low value for 
the radio flares. Because a single-zone model only contains a sin¬ 
gle Doppler factor, this parameter would need to change rapidly 
between the two flares. This scenario would require the Doppler 
factor to drop by > 10 between the y-ray and radio flare within a 
period of 4 = 60 d, which is much larger than the modest variations 
applied in the modeling of the 2012 flares. 

The 2012 y-ray flare had two prominent spikes, and we used 
these to investigate the physical conditions required to produce 
the extreme radio event. By using the 2013 model parameters as 
a starting point and introducing as few modifications as possible, 
we considered two possible scenarios. These result from the as¬ 
sociation of the radio flare with the first or second y-ray spikes. 
Under the assumption of a physical connection between the radio 
flare and the first y-ray spike, we showed that the addition of a 
varying Doppler beaming factor to a one-zone SSC model with 
fairly typical parameters can explain the observed sharp radio flare. 
The SSC model would otherwise produce too broad a radio flare, 
similar to the 2013 event (see Figs. I4l5l >. Although this scenario 
succeeds in explaining the radio flare, it results in a wider y-ray 
pulse profile than is observed (solid lines in Fig. O. Moreover, it 
is difficult to explain the lack of a second peak in the radio light 
curve, as discussed in the end of Sect. 4.2.1. Nonetheless, a varying 
Doppler beaming factor is not unreasonable. An adequate varia¬ 
tion could result from, e.g., a modest change in the Lorentz factor 
or the viewing angle. In the latter case, the sharp flare would oc¬ 
cur when the direction of motion of the emission region crosses 
very near to the line of sight to the observer. Similar models with 


blazars as well (e. 

e.. Marscher et al.l2008l. 2010l; Abdo et alj2010al; 

iRani et al.l l2013l; 

Aleksic et alj 2014; Molina etalj 2014). These 


models assume that either the emission feature mo ves on a helical 
trajectory due to an ordered hel ical magnetic field jMarscher et al.l 
l 2008 l I2OIOI: iMolitra et al ' 

the jet l lAbdo et al .II2OIO1 


I] |2014h. or that there i 

^: lAleksic et al ][20li) . 


is an actual bend in 


In our second scenario, we considered the possibility that the 
radio flare is associated with the second, broader y-ray spike. In 
this case, the short duration of the radio flare and the shorter time 
delay between the radio and y-ray flares imply faster electron cool¬ 
ing, and thus a stronger magnetic field than the one used in the first 
scenario. We found indeed that a viable model is achieved by in¬ 
creasing the magnetic field strength by a factor of 2.5. Fig.|5]shows 
that the model light curve (dashed lines) describes well the y-ray 
data, although the modeled radio flare is now not quite as sharp 
as the observed one. In this scenario, the first y-ray spike would 
have to be unassociated with the radio event, perhaps by occurring 
closer to the black hole where the emission region is optically thick 
for radio emission. 


6 CONCLUSIONS 

We have obtained radio, millimeter and y-ray light curves of 
Mrk 421 during the flaring activity in 2012 and 2013. In July 2012, 
Mrk 421 exhibited the largest y-ray flare observed by Fermi since 
the beginning of its mission. About 40 - 70 d later, the largest ever 
15 GHz flare was observed. The flare rise time determined from an 
exponential fit was just 10.6 ± 0.5 d, which is extreme compared 


to previous radio flares observed in the source. This implies a vari¬ 
ability Doppler factor 6 ~ 3 - 10. In 2013 Mrk 421 underwent 
major y-ray flaring again, followed by radio and millimeter flares 
about 60 d later. Under the assumption that the events in the differ¬ 
ent bands are physically connected we have modeled the variations 
with the simplest possible theoretical model, with the main goal of 
explaining the extreme radio flare in 2012. Starting with the less 
extreme 2013 flare, we obtained a one-zone SSC model that could 
explain the main features of the flares reasonably well. The 2012 
y-ray flare was double peaked, and by modeling the radio flare as 
connected to each peak in turn, we were able to reproduce the radio 
flare with either a varying Doppler factor, or an increased magnetic 
field strength. In both cases several specific conditions in the jet 
need to be fulfilled for the models to be viable, showing the ex¬ 
treme and unique nature of the 2012 event. 
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